Huge anisotropic magneto-resistance in iridium atomic chains 
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We analyze in this article the magneto-resistance ratio of finite and infinite iridium and platinum chains. 
Our calculations, that are based on a combination of non equilibrium Green function techniques and density 
functional theory, include a fully self-consistent treatment of non-collinear magnetism and of the spin-orbit in- 
teraction. They indicate that, in addition to having an extremely large magnetic anisotropy that may overcome 
the super-paramagnetic limit, infinite and also realistic finite-length iridium chains show sizeable anisotropic 
magnetoresistance ratios. We therefore propose iridium nanostructures as promising candidates for nanospin- 
tronics logic devices. 
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The ability to enhance and tailor at the same time the mag- 
netic anisotropy and magneto-resistance of atomic-sized mag- 
netic bits and junctions will determine whether nanospintron- 
ics will be a viable technology. Cobalt clusters and chains 
with magnetic anisotropy barriers of the order of 1 meV were 
fabricated a few years ago by deposition on substrates of 5d 
elements iH 01 ■ Furthermore, several recent theoretical predic- 
tions have pointed out that atomic clusters and chains made 
of 4d elements should have even higher magnetic anisotropy 
barriers than Co nanostructures Q 0. Bt). In the case of 5d 
elements like platinum and iridium, calculations predict that 
the anisotropics of iridium and platinum atomic clusters and 
chains can be even larger than room-temperature thermal en- 
ergies, so that these nanostructures can overcome the techno- 
logically relevant superparamagnetic limitfl S 01 ■ 

Suspended chains of gold atoms connecting two gold elec- 
trodes were fabricated using the scanning tunneling micro- 
scope and mechanically controllable break junction (MCBJ) 
techniques^ where a quantized conductance close to 
Go = 2 e^//i was measured. Subsequent experiments found 
that platinum and iridium also form suspended chains, but that 
other elements like Ni, Co, Pd or Rh do not un- 
less doped with oxygen or other elements jlU. It was also 
demonstrated that longer atomic chains could be made by de- 
position on stepped surfaces lit], or by encapsulation in carbon 
nanotubesll3l. 



The reduced dimensionality of nanometric objects in- 
creases their tendency towards magnetism. Ugarte and 
coworkers proposed that their experimental conductance 
data for platinum chains could be explained provided that 
these were magnetic 11 1411 . Subsequent theoretical work in- 
deed found that 4d and 5d infinite linear chains of bulk- 
paramagnetic materials do become magnetic 1 15 Till . Very re- 
cently, several groups have found that the magnetic anisotropy 
energy (MAE) per atom of infinite zigzag platinum and irid- 
ium chains j6[l, and of linear platinum chains f?! show values 
as large as 50 meV to 100 meV when stretched, which cor- 
respond to equivalent temperatures of 500K to 1000 K. This 
suggests that the blocking temperature Tb of such nanostruc- 



tures will attain a similar value and therefore they should over- 
come the superparamagnetic limit, so that their magnetization 
would be finite at room temperature. 

We draw the attention in this article to iridium and platinum 
atomic structures as promising candidates for nanospintron- 
ics. We show that iridium atomic chains not only have huge 
magnetic anisotropics, but also meet the requirement of hav- 
ing a significant magneto-resistance ratio. The simulations 
presented here are based on density functional theorv 11711 
as implemented in the SIESTA code llSll . and include a 
fully self-consistent and spin non-collinear implementation 
of the spin-orbit interaction lfigll . This spin-orbit implemen- 
tation has been successfully used to compute the magnetic 
anisotropy of infinite Ir and Pt atomic chains and of tran- 
sition metal atomic clusters 1411, in agreement with the re- 
sults of other simulations Isllsl 01 • The conductance cal- 
culations were carried out with our spin non-collinear code 
SMEAGOL|201, which includes the fully self-consistent im- 
plementation of the spin-orbit interaction referred above. We 
note that SMEAGOL has been shown to provide conductance 
data that compare rather accurately with the available experi- 
mental data on platinum chains lid lul |2l|l. The simulations 



have used the Local Density Approximation |221 and norm- 
conserving pseudopotentials 12311 . The atomic parameters that 
were fed in the pseudopotential generator, as well as the basis 
set used are similar to those reported in references S 21 1 ■ 
Figure 1 (c) shows a schematic view of the experimental 
setup of a MCBJ experiment, where at the last stages before 
breakage of a Pt or Ir strip, a short suspended atomic chain 
is formed. In these experiments, an electric current is applied 
to the circuit, and the conductance G is measured many times 
as the chain elongates. The whole setup can also be under- 
stood as a magnetic junction that connects two paramagnetic 
electrodes. Recent work has shown that theory can reproduce 
accurately the experimental conductance data of finite-length 
platinum chains 12111 . but has not found significant differences 



between the paramagnetic and the magnetic cases 112411 . 

We show the energy of linear platinum and iridium chains 
as a function of in Figs. 2(a) and 3(a). It is important to 
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FIG. 1 : (Color online) The geometry of the chains simulated in this 
article, referred to the three coordinate axes, (a) Ladder geometry 
that corresponds to the absolute minimum of the infinite chains, see 
the top right panel in Figures 2 and 3; (b) zigzag geometry that cor- 
responds to the second minimum of the infinite chains, see the top 
right panels in Figures 2 and 3; (c) schematic view of the geometry 
of the relaxed finite-length chains in contact to two electrodes. 



notice that these chains are metastable. Indeed, allowing for 
lateral displacements of the atoms yields the energy curves 
in Figs. 2(b) and 3(b), that only have minima for ladder or 
zigzag arrangements similar to those shown in Figs. 1(a) and 
(b). Furthermore, relaxation of the forces always finishes in 
zigzag or ladder chains, unless the atoms are constrained to 
lie in a linear array. Actually, ladders are energetically more 
stable than zigzag geometries. But the relaxation of the forces 
in a MCBJ geometry, like the one shown in Fig. 1(c) only 
leads to zigzag arrangements, that gradually straighten as the 
electrodes are pulled apart. Related to this effect, we note that 
infinite platinum (iridium) zigzag chains become linear only 
for distances longer than 2.55 A (2.40 A), since at these 
distances, the zigzag angles fall abruptly to zero|0]. These 
values of are marked by a dotted vertical line in the right 
panels of Figures 2 and 3. We will therefore only present re- 
sults for the conductance of infinite zigzag and linear, as well 
as finite zigzag, chains in this article. 

We have performed separate simulations for orientations 
of the atomic spins along the chain axis (z- or easy-axis), 
perpendicular to it, but still in the zigzag plane (x-axis) and 
perpendicular to the zigzag plane (y-axis), but the energies 
Ex,y,z can not be discriminated in the energy scale used in 
Figs. 2(a,b) and 3(a,b). We have therefore computed!^ the 
energy differences A^^^j, = — E^.y per atom, that we de- 
fine as the MAEs, and plotted them in Figs. 2(c,d) and 3(c,d) 
for platinum and iridium respectively. Our definition implies 



FIG. 2: (color online) Results of the simulations of infinite platinum 
chains as a function of the distance dz defined in Figure 1 for orienta- 
tions of the spins along the = x,y and z axes. The left panels corre- 
spond to the case of metastable linear chains. The right panels corre- 
sponds to the geometries (b) and (c) in Figure 1 . The dashed vertical 
line indicates the distance dz at which the zigzag chains cease to be 
stable, (a) Energies E^.y.z per atom; (b) magnetic anisotropy ener- 
gies per atom Ax,y', (c) spin moments per atom; (d) conductance. 



that when both A^. j, are negative, the magnetization lies par- 
allel to the easy axis. On the contrary, when Ax,y is positive, 
the magnetization is oriented perpendicular to the easy axis, 
and aligns either parallel to the x- or the j/-axes, depending 
on whether A^; is larger or smaller than Aj,, respectively. No- 
tice that our results for linear platinum chains (Fig. 2(c)) are 
similar to those presented in Ref. 101. More interestingly, the 
anisotropy of zigzag platinum and iridium chains is x-axis for 
short dz, but shifts to easy axis if the chains are elongated, 
with the MAE achieving values as large as 50-80 meV de- 
pending on the specific value of dz- The third row in Figs. 
2 and 3 shows the spin moment per atom as a function of dz- 
Smogunov and coworkers||3l found that linear platinum chains 
become magnetized at a shorter elongation dz for spin lying in 
the axis of the chain than for spins oriented perpendicular to 
it. This fact led them to predict a colossal magnetic anisotropy 
for the window of elongations within the onset of both magne- 
tizations. However, we note here that the magnetization of the 
more realistic zigzag chains is finite for almost the full range 
of elongations regardless of the spin orientation, and therefore 
that effect is not expected to be seen. Interestingly, we find 
that the magnetization of iridium chains decreases when the 
chain is pulled apart, in a window of elongations around the 
equilibrium dz ■ 

Since infinite Ir and Pt chains show large MAEs, an obvi- 
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FIG. 3: (Color online) Results of the simulations of infinite iridium 
chains as a function of the distance dz defined in Figure 1 for orienta- 
tions of the spins along the = x,y and z axes. The left panels corre- 
spond to the case of metastable linear chains. The right panels corre- 
sponds to the geometries (b) and (c) in Figure 1 . The dashed vertical 
line indicates the distance at which the zigzag chains cease to be 
stable. )a) Energies E^.y.z per atom; (b) magnetic anisotropy ener- 
gies per atom Ax,y', (c) spin moments per atom; (d) conductance. 

ous question is whether they also show a strong anisotropic 
magneto-resistance effect that could be activated by the pres- 
ence of a magnetic field. If so, this effect could open the way 
towards room-temperature inorganic nanomagnetic memory 
bits. We have computed the low- voltage conductance Gx,y,z 
for orientations of the spins along the three axes, that we show 
in Figs. 2(g,h) and 3(g,h), again as a function of dz- We have 
found that the magneto-resistance ratio, defined as 

AMR,,y = 2 (1) 

is finite but small over wide ranges of elongations dz- Further- 
more, iridium chains show larger magneto-resistance ratios 
than platinum ones for a wider range of dz at and above the 
equilibrium elongation. We therefore expect that the magneto- 
resistance of iridium chains should be more robust than that of 
platinum chains. 

Remarkably, for certain windows of the elongations dz in 
zigzag chains, one of the conductances may fall exactly to 
zero, marking the onset of a perfect magneto-resistance ra- 
tio. Notice though that these windows of perfect magneto- 
resistivity do not appear for linear chains. Zigzag platinum 
chains show three of those windows, which occur when the 
atomic spin moments are oriented in the y-axis (i.e. perpen- 
dicular to the zigzag plane). Two of these appear for unreal- 
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FIG. 4: Band structure of a zigzag platinum chain with an elongation 
dz = 2.45 A. From left to right the panels show the bands for a 
paramagnetic calculation, and for spin-orbit simulations whereby the 
spins are oriented along the z-, x- and y-axes, respectively. 



istically short elongations d^, but the third develops for dis- 
tances slightly larger that the equilibrium d^ and therefore 
could have the potential to be realized and measured exper- 
imentally. Iridium chains show only one of these windows, 
where Gx vanishes, which occurs for elongations about the 
equilibrium dz , where the magnetic anisotropy favors orienta- 
tions of the spins along the x-axis. The MAE in this window 
of elongations is of the order of 25 me V. Interestingly, the van- 
ishing conductance develops precisely for orientations of the 
spins along the x-axis. 

Since the above finite magneto-resistance ratios originate 
from the different number of channels available for electron 
conduction close to the Fermi energy for the different spin ori- 
entations, we now study the band structure of infinite chains. 
Indeed, the number of conductance channels can be found by 
counting the number of bands crossing the Fermi energy. We 
note that each band moves up or down in energy at a different 
pace as the chain elongates. Furthermore, since the spin-orbit 
interaction couples spin and orbital degrees of freedom, these 
shifts are different for the different spin orientations. Since the 
band structure of non-conducting chains must show a mini- 
gap at the Fermi energy, we have plotted in Fig. 4 the band 
structure of a platinum chain at dz — 2.45 A for paramag- 
netic as well as for spin-orbit simulations. We indeed find that 
when the atomic spins point in the y-axis there exists a small 
mini-gap Eg = 25 meV atk = OAn/dz- This value allows us 
to estimate the minimum number of atoms that a finite-sized 
platinum chain must have to develop such a mini-gaps. Us- 
ing the relationship ~ hvp/Eg dz, where we estimate the 
Fermi velocity at this fc-point ashvp — 1000 meV dz , we find 
a critical length of about 40 atoms. We therefore expect that 
it will be difficult to measure a significant magneto-resistance 
ratio in platinum chains. We have indeed simulated 5-atoms 
long platinum chains contacted to (001) platinum electrodes, 
in a geometry similar to Fig. 1 (c) and found that the magneto- 
resistance ratio was essentially zero. 

In contrast. Fig. 5 shows the band structures of iridium 
chains, obtained from paramagnetic and also for spin-orbit 
simulations for an elongation dz such that Gx falls to zero. We 
note that the spin-orbit interaction generates much larger gaps 
for iridium than for platinum. For this particular elongation. 
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FIG. 5: Band structure of a zigzag iridium chain with an elongation 
dz = 1.96 A. From left to right the panels show the bands for a 
paramagnetic calculation, and for spin-orbit simulations whereby the 
spins are oriented along the z-, x- and y-axes, respectively. 
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FIG. 6: (Color online) Conductance as a function of of iridium 
chains containing 21 atoms. 

the gap when the spins point along the x-axis is positioned 
exactly at the Fermi energy and has a value Eg ~ 100 meV at 
the r point. Using the estimated Fermi velocity hvp — 1500 
meV dz, we predict that chains containing about 15 atoms 
should show large magneto-resistance ratios. 

To confirm this, we have simulated 2- and 21 -atoms long 
iridium chains in contact to (111) iridium electrodes termi- 
nated by a pyramid, as in Fig. 1 (c), for a number of elonga- 
tions of the chain, which cover both the equilibrium distance 
and also stretched situations. We have found that the magneto- 
resistive ratio of the 2-atom chains was negligible for whatever 
of these elongation. In contrast, we have found that Gx is in- 
deed substantially smaller than Gy^z, as we show in Fig. 6, 
rendering a finite ratio AMRx- We stress that our previous 
analysis of Figs. 3 and 5 indeed led us to expect that Gx had 
to be smaller than Gy,z. 

To summarize, we expect that short platinum chains should 



show a negligible magneto-resistance ratio until the number 
of atoms in the chain exceeds of order 40. Since such a 
length is very difficult to attain experimentally, we believe 
that platinum is not the optimum magneto-resistance material. 
We predict that a better candidate element is iridium. Short 
iridium chains are expected to display small but measurable 
magneto-resistance ratios. Longer chains could develop per- 
fect magneto-resistance for certain elongations. These clear 
differences are experimentally accessible and therefore we an- 
ticipate that future experimental studies of iridium chains will 
prove to be extremely fruitful. 
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